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Abstract 

The level dynamics of pseudointegrable systems with different genus numbers g is 
studied experimentally using microwave cavities. For higher energies the distribu- 
tion of the eigenvalue velocities is Gaussian, as it is expected for chaotic systems 
with time-reversal symmetry, and shows no dependence on g. Also the curvature 
distribution P(k) for large k is decaying as it is expected for chaotic systems, i.e. 
P(k) ~ |&| -3 - For small k an intermediate behavior is found, where P(k) changes 
from integrable towards chaotic behavior with growing g. 
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1 Introduction 



The motion of a classical particle in a billiard system can show regular, chaotic 
or intermediate behavior, depending on the billiard geometry. A potential 
well of the same geometry as the corresponding classical billiard - a quantum 
billiard - reflects this regular, chaotic or intermediate behavior in the statistics 
and the dynamics of its eigenvalues. From spectra only static correlations like 
nearest neighbor distributions can be obtained [1] . Level dynamics additionally 
observes eigenvalues in dependence of an external parameter. 
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Pechukas [2] and Yukawa [3] established the formal analogy between a one- 
dimensional classical gas and a perturbed quantum system H = H +tV, where 
Hq is the Hamiltonian of the unperturbed system, and tV is a perturbation. 
Here, the change of the eigenvalues E n (t) with the parameter t is governed 
by the same "equations of motion" as the positions x n (t) of the gas particles 
propagating with time t. In this picture, a chaotic system can be described 
as a one-dimensional gas of particles with repulsive long-range interactions 
between them. An integrable system, on the other hand, refers to an ideal gas 
with no interaction or one close to zero, when tunneling effects are taken into 
account [4]. The same equations of motion as for the Pechukas- Yukawa gas 
have been derived for the billiard system under the shift Al of a billiard wall 
[5,6], where Al (or the length of the billiard /) refers to the time t. 

In this paper, we are interested in the level dynamics of two-dimensional pseu- 
dointegrable billiards [7]. Examples are polygons with only rational angles 
riiir/mi, with n^rrii G N and at least one n» > 1. They are characterized by 
their genus number g — 1 + (M/2) J2'i=i( n i ~~ where J is the number 

of angles and M is the least common multiple of the rrij. The level statistics 
of pseudointegrable systems has been found intermediate between those of 
integrable and chaotic systems [8,9,10,11,12,13,14]. In the following, we con- 
centrate on their level dynamics, i.e. on the distribution of the velocities and 
curvatures of the eigenvalues E n (l). 

It has been derived in [3] that the distribution of the velocities v n = dE n /dl 



of chaotic systems is Gaussian, P(v) = y ^ exp (— |t> 2 ). The level curvature 
is defined as the second derivative of the eigenvalues K n (l) = d 2 E n (l)/dl 2 . We 
scale it to the dimensionless curvature 

KA 

<K<{dE n /dl) 2 >' (1) 



where A denotes the average level spacing [15,16,17]. The trajectories of parti- 
cles without interaction are straight lines, i.e., all k n (l) = 0, which corresponds 
to P(k) ~ 5(0) for integrable systems. In a chaotic billiard, on the other hand, 
the interaction between the "eigenvalue particles" are long-ranged, leading 
to curved trajectories with a broad distribution P(k). The curvatures k for 
pure random Gaussian orthogonal ensembles (GOE) [1] obey the generalized 
Cauchy distribution [15,16,17] 



Eq.(2) corresponds to a power-law behavior for large k, P(k) ~ \k\ 3 , in 
agreement with earlier derivations [18]. The small curvature behavior has been 
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Fig. 1. Investigated billiards with genus g = 3 and g = 12. One wall was shifted in 
steps of 0.25 mm. Some of the widths w^bo and the length i^bo of the main bouncing 
ball orbits, taken into account for the unfolding of the spectra, are shown. 

found to be non generic, when families of periodic orbits like bouncing balls 
are not taken explicitely into account [15,19] (see below). 

Numerical works on pseudointegrable billiards, as the barrier billiard [20] and 
systems approaching the chaotic Sinai billiard by polygons with different genus 
numbers between g — 2 and g — 7 [21] found that the trajectories of the eigen- 
values do not cross each other (as in integrable billiards), but show "avoided 
crossings", indicating that some repulsive interaction must be present. In [22], 
a model of short-range interactions was used to calculate the spectral statistics 
of pseudointegrable systems. In [21], the curvature distribution P(K) was cal- 
culated, which showed the expected GOE behavior for large K. Supprisingly, 
a peak was found for \K\ close to zero that was larger in the case of the Sinai 
billiard than for the polygon billiards. As also pointed out by the authors, this 
demands further investigation and we will show that it might be due to the 
fact that the bouncing ball orbits between parallel walls were not taken into 
account explicitely. 



2 Systems and Experiment 

In this work, we measure experimentally the eigenvalues E n of pseudointe- 
grable microwave billiards under the shift of a wall [23]. We investigate two 
pseudointegrable systems with genus g = 3 and 12 (see Fig. 1). The widths 
of both systems are equal 237.1mm and the initial lenghts are Iq = 380 mm. 
The length of the system / can be changed by shifting one wall, which is used 
as the external parameter. The wall was moved in steps of 0.25 mm. For each 
position of the wall a spectrum from 0.5 GHz up to 18.1 GHz was measured. 
In total we have 161 measurements giving 4cm as a total change of length. A 
detailed analysis of all obtained spectra was performed up to a frequency of 
15 GHz, yielding about 400 resonances for each billiard. For higher frequencies 
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Fig. 2. The unfolded energy levels i = 151, . . . , 180 as a function of the billiard 
length I for the g = 3 -system (a) and the <? = 12 -system (b). It can be recognized 
that there are more straight lines in the spectrum of the billiard with g = 3. 

the resonances could not be distinguished with sufficient reliability. 

In order to perform a statistical analysis of the given eigenvalue sequences inde- 
pendently from the change of size, the measured spectrum has to be unfolded 
to a mean level spacing of 1. To this end we used the Weyl formula [24] with 
the additional correction terms N bbo (E) that bring into account the existence 
of bouncing ball orbits (bbo) between parallel segments of the system, 

e n = N(E n ) = ^-E n - ^-^ l + Y / N bbo (E n ), (3) 

i 

where A is the area and B the perimeter of the billiard. The correction term 
derived in [25] 

N- bo (E n ) = ^%\[^E^ £ cos(2ml bbo , t JK - ^), (4) 

(2tt) ' V Zbbo,i m= i 4 

is used additively for every parallel segment in our systems. Here Wbbo,i is 
the width of the i-th parallel part of the billiard boundary and £bbo,i is the 
distance between its opposite walls (see also Fig. 1). From now on, we refer to 
the unfolded energy values e n . 



3 Results 

In Fig. 2 we present a part of the unfolded energy spectra e n as a function of 
the length / of the billiard for both systems. Both figures look quite similar, 
though the amount of straight lines seems to be higher for the g = 3 -system 
than for the g — 12 -system, which could be due to a shorter range of the 
interactions between particles in pseudointegrable systems. In the following, 
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Fig. 3. Distribution P(v) of the rescaled velocities v for systems with g = 3 and 
g = 12 in the high energy region (e G [350,400]). The solid line is a Gaussian. 

we investigate this assumption quantitatively by calculating the velocity and 
the curvature distribution functions. 

In the case of GOE systems [2,3], one expects a Gaussian distribution for the 
velocity distribution P(v) of the scaled velocity v = (de/dl)/ < (de/dl) 2 > 1 ^ 2 , 
P(v) = (l/-\/27r) exp (— v 2 /2). However, for low-lying energies, the statisti- 
cal properties often deviate from theoretical expectations. Therefore we in- 
vestigated P(v) for the low-energy and the high-energy regime separately. 
Fig. 3 (a, b) shows the high-energy regime, e G [350,400], which is for both 
systems close to the Gaussian curves of chaotic systems (whereas large devi- 
ations occur at low energies). 

For studying the curvatures k we chose the same energy interval e G [350, 400] 
as in Fig. 3. Fig. 4 (a, b) shows P(k) for large k in a double-logarithmic scale 
for g = 3 and g = 12. Both experimental results show a |/c|~ 3 behavior of the 
tail (see solid line), which is expected for chaotic systems (GOE). Therefore, 
it seems that the decay is not depending on the genus number. Similar results 
were already reported for the pseudointegrable billiards of [21]. This can be 
understood as follows: The large curvature behavior is mainly governed by 
the curvatures at the close encounters i.e. the avoided crossings. There, the 
curvature only depends on the interaction of the eigenvalue with its nearest 
neighbor. At short distances, the exact range of the interaction is unimportant 
and already short-range interactions are sufficient to lead to the same |A:|~ 3 
behavior as in GOE systems. 

The distributions of small values of k, on the other hand, may depend on 
the range of the interactions. In Fig. 4 (c), we therefore show P(k) on a 
linear scale and compare it to the theoretical prediction for a GOE system 
(solid line, cf. Eq. (2)). Whereas for the g = 3-system, the values of P(k) 
at small k are much higher, P(k) of the g = 12-billiard is indistinguishable 
from Eq. (2). This corresponds to the observation of more straight parts of 
eigenvalue trajectories in the g = 3-system (see Fig. 2) and indicates that 
P(k) changes from a P(k) ~ 5(0) for integrable systems towards Eq. (2) with 
growing g. 
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Fig. 4. Distribution of the absolute value of the scaled curvatures |fe| in a double 
logarithmic scale for the g = 3 -billiard (a) and the g = 12 -billiard (b). The solid 
line indicates a |&|~ 3 behavior as expected for GOE systems. Case (c) shows the 
distribution of k for g = 3 (circles) and g = 12 (triangles) for e £ [350,400]. The 
solid line shows the theoretical prediction for a GOE system (see Eq. (2)). 



4 Summary and Conclusions 



In summary, we investigated experimentally the level dynamics of two pseu- 
dointegrable billiards with different genus numbers g and calculated the dis- 
tribution P(v) of the velocities and P(k) of the curvatures of the eigenvalues. 
While the velocity distribution as well as P(k) for large k shows the behavior 
expected for GOE systems, P(k) depends on g. Using an unfolding procedure 
that takes the contributions of the bouncing ball orbits between parallel seg- 
ments of the system into account, we found that P(k) approaches Eq. (2) of 
the GOE system with increasing g. It will be very interesting to investigate 
also systems with 3 < g < 12 as well as with other system details, to see if 
this behavior is general. 

The distribution P(k) for small curvatures does not correspond to the results 
presented in [21], which could be due to the different unfolding procedure. 
Contrary to our work, the unfolding by [21] did not take into account the 
bouncing ball orbits between parallel walls and it is not evident what governs 
the behavior of the small curvatures in this case. 

Our results (see Fig. 4 (c)) support the model of short-range interactions 
between eigenvalue particles in pseudointegrable systems. Small curvatures 
arise from eigenvalues in large distances from each other. When the interaction 
between them decays rapidly, eigenvalues beyond a certain distance do not 
influence each other and their trajectories become nearly straight lines. Our 
results indicate that the range of the interaction increases with g. 

We thank H.-J. Stockmann and R. Schafer for valuable discussions and the 
Deutsche Forschungsgemeinschaft for financial support. 
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